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Abstract. We study numerically the effects of static imperfections and residual couplings between qubits for the 
quantum phase estimation algorithm with two qubits. We show that the success probability of the algorithm is af- 
fected significantly more by static imperfections than by random noise errors in quantum gates. An improvement of the 
algorithm accuracy can be reached by application of the Pauli-random-error-correction method (PAREC). 

PACS. 03.67.Lx Quantum computation - 85.25. Cp Josephson devices - 24. 10.Cn Many-body theory 



1 Introduction 

Quantum computers [1| are doing steady progress increasing 
the number of qubits and accuracy of quantum gates. Among 
most advanced physical implementations with possible scal- 
able architecture are ion based quantum computers (see e.g. j2| 
and Refs. therein) and solid state superconducting qubits (see 
e.g. EJUS]). In the present situation when only a few qubits 
are available and the quantum gate accuracy is limited it is in- 
teresting to test a performance of simple quantum algorithms 
operating at such realistic conditions. One of such algorithms 
is the Iterative Quantum Phase Estimation Algorithm (IQPEA) 
proposed recently by Dobsfcek et al. 0. According to the re- 
sults obtained there the IQPEA works reliably even in presence 
of relatively strong noise in quantum gates. The algorithm [6 1 is 
based on the semiclassical Quantum Fourier Transform (QFT) 
which uses one ancilla qubit, iterative measurements and 
a classical information feedback. The advantages of the algo- 
rithm are the following: it uses only a single ancillary qubit and 
its theoretical accuracy of the eigenvalue found is limited only 
by the number of times the algorithm is applied. Due to that the 
IQPEA can be used as a benchmark algorithm for the maximal 
accuracy that could be obtained in a given experimental setup 

We describe briefly the IQPEA in the case of a two qubit 
circuit proposed in [6 1 as a minimal benchmarking circuit. The 
goal is to measure the eigenphase <j> of some operator U with 
precision set to m significant bits. The standard Phase Estima- 
tion Algorithm (PEA) [8 1 requires m ancillary qubits, to get the 
desired precision, and the possibility to implement efficiently 
control-f/ 2 gates. Using the semiclassical implementation of 
the QFT an alternative algorithm with only one ancillary qubit 
and measurements can be designed |9|. One further step ahead 
is done in [6| where a feedback of the measurement result is 
used to correct the phase. In this way the IQPEA provides a 
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way to compute the phase theoretically with arbitrary preci- 
sion. The proposed scheme can be used as a benchmarking cir- 
cuit which is tolerant to a rather strong random noise in quan- 
tum gates. 

In 1 6 1 only the case of random noise errors in quantum gates 
and environment dephasing are considered. At the same time it 
is known that a presence of static imperfections and residual 
couplings between qubits may lead to an emergence of quan- 
tum chaos in a quantum computer hardware iflOl . Such static 
imperfections affect the accuracy of quantum computation in a 
significantly stronger way compared to random errors in quan- 
tum gates 11 1 Ul 21 . Thus it is interesting to test the effects of 
static imperfections in IQPEA with a small number of qubits, 
e.g. two qubits. Indeed, our studies presented in this paper show 
that the static imperfections lead to a significant drop of the 
computation accuracy and the algorithm success probability. To 
correct these quantum errors induced by static imperfections in 
IQPEA we apply the Pauli-random-error-correction (PAREC) 
method proposed in [ 1 3 ] and tested in various quantum circuits 
B14II15L Our results for IQPEA show that the PAREC allows to 
improve significantly the accuracy of quantum computation. 

The paper is organized as follows. First we briefly describe 
the IQPEA for a 2-qubit system (Section 2). We then compare 
the effects of random phase errors in quantum gates and the 
effects of static imperfections (Section 3). Then we show how 
the PAREC method corrects the errors induced by static im- 
perfections (Section 4). The summary of the results is given in 
Section 5. 



2 Brief description of IQPEA 



The goal of IQPEA is to find the eigenphase of an operator U. 
We consider the simplest case with one-qubit unitary diagonal 
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operator in the computational basis 

U- 



— i2?T0 q 

e i2 ^ 



(1) 



with (j) e [0, 1]. We want to find (j) with up to m bits of accuracy 
(or error smaller than 2~"'). The IQPEA |6| procedure consists 
in applying m times the circuit shown in Fig. Q]to the 2-qubit 
state 1 00). For simplicity we use the presentation 



H f 

= 5>-2-'' = O.0ifc...fc n OOO. 



(2) 



assuming that the binary expansion of is finite. For the first 
step (i = 0) we take o, = so the Z-rotation does not act. After 
this first run 



1 00 



1 



? i2 ^)|00) + (l-e l2 ^)|10) 



(3) 



and the measurement of the "left" qubit yields Po(|0)) = 
cos 2 (7r(O.0 m )) which is unity if <j) m = and zero if 0,„ = 1. Thus 
the least significant bit of <p is obtained deterministically. The 
key element is that in the following steps of the algorithm we 
use the classical information obtained from the measurement 
to correct the phase by a Z-rotation. Before the last Hadamard 
gate the phase in the second step is 2n(0.(p m - 1 0,„OO . . .) and af- 
ter performing a Z-rotation with Ofy = — 27r(O.O0 m _ i ) the prob- 
ability becomes Pi(\0)) — cos 2 (7r(O.0 m _i)). Consequently, the 
result of the first measurement is used as a feedback for the 
algorithm to obtain the second least significant bit which is 
obtained deterministically. In theory, following this procedure 
each bit can be obtained. 



|0> -H 
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Fig. 1. Step (where i = 0, . . . , m- 1) of the IQPEA of (6). Here R z 
is a rotation by an angle ft), around z axis. 



In reality the phase is 

<j) = + 82' 



(4) 



where = 0.0102 • • • m OOO and 8 E [0, 1) is the reminder. The 
probability of measuring 0,„ correctly is thus 

P 1 (|0»=cos 2 W(0.^) + 5/2))l 2 
P 1 (|l))=sin 2 (rc((0.<M + S/2)) J- Pl - C0S W 2 )- (5) 

The next step gives P2 = cos 2 (7r5/2 2 ) and eventually P^ = 
cos 2 (n 8 /2 k ) so that the total probability of measuring the phase 
correctly is given by 



P tot (S) = Y\cos 2 (nS/2 k ) = 



sin 2 (7r5) 



k=l 



2 2m sm 2 (n8/2 m ) 



(6) 



The success probability in (O is bounded in the limit m — * °° by 
4/7r 2 [8|. In fact, the rounding error permits us to neglect the 
least significant bit and consider as probability of success the 
sum Ptot (<5)+P tot (l -8) = 8/tt 2 when yyi ► 00 1 6 1 . This lower 
bound could be raised by repeated measurement of the first few 
bits and majority vote [8]|6). 



3 IQPEA and static imperfections 

We consider two kinds of circuit imperfections: random phase 
errors in rotations and static imperfections due to residual cou- 
plings between qubits. To model random quantum phase errors 
we assume that the rotation on angle 9 



R aiv) (9)=exp[-ia( v) d/2 



(7) 



(with a( v ' a Pauli operator) is replaced by rotation on angle 
0(1 + A ) with A randomly and uniformly distributed in the 
interval 

^[-|,|]. CD 

In other words the original rotation Hamiltonian has now an 
additional term 

8EU = ^oM (9) 

Each gate in Fig.[T]is implemented with rotations having differ- 
ent random realizations of A . This is the case of random noise 
errors considered in [6 | where it was shown that the algorithm 
is rather robust. 




e i- £ 2 



Fig. 2. (Color online) Success probability for the algorithm to de- 
termine the phase with a precision of up to 10 bits, as a function of 
the parameter e = £j = £2 characterizing the error strength. Symbols 
mark: random phase errors in the Hadamard gates(D), the R z gate (■), 
the controlled-C 2 (o) and errors in all of the gates (•). The case of 
static imperfections is shown by (A)/(A) in absence/presence of ran- 
dom phase errors in the gates. Averaging is done over 2000 randomly 
chosen phases. 



On the other hand, we model the effects of residual static 
couplings between qubits by an imperfection Hamiltonian of 
the form used in 01011121 : 



8H M (x) = 8ig\ z> + 820™ +2Jaf'a : 



(x) 



(10) 
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where a- are the Pauli operators acting on the iih qubit and 
8(, J are random coefficients uniformly distributed according to 

Si, J G [— a V3e 2 ,aV3e 2 ] , (11) 

(with a constant). We suppose that between each gate in the 
algorithm there is a finite time At which remains fixed during 
the algorithm and that 8H acts via the unitary propagator 

U ihll = e' 5H (12) 

where the time At has been absorbed into the constants 8 and 
J in ([Tol l. The time At can be considered as an effective gate 
duration, a similar scheme is used in lfl2l . 

In order to compare the effects of both types of errors we 
compute (tr[8H 2 ]) e 2 . The value a « 0.37 is determined so 
that 

(tx[SHi}) « (tx[SHl}) if e, = e 2 , (13) 

the approximation is done taking 9 = % in in Eq. (0. 

The static type of imperfections is especially important since 
generally the errors produced in this case are accumulated co- 
herently that leads to a quadratic term in the decay of fidelity 
1 12 1 thus limiting considerably the maximum time over which 
an accurate quantum computation can be performed. 

In Fig.|2]we show the success probability of measuring cor- 
rectly the phase 0, for a chosen accuracy of 2~ 10 , as a function 
of the parameters E\ = £ 2 = £■ We averaged over 2000 uni- 
formly random phases (j) £ [0, 1]. For the cases where only ran- 
dom phase errors act we see that the algorithm is rather robust. 
Indeed, the decay of the success probability is relatively slow 
even when all the gates involved have errors. On a contrary, a 
dramatic drop of the accuracy of computation is seen when we 
include the effects of static imperfections. The decay is much 
faster than for random phase errors even for comparable values 
of tr[8H 2 } corresponding to a typical experimental situation. 
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Fig. 3. (Color online) Schematic representation of the PAREC method. 
The vertical thick (red) lines indicate the place where the static cou- 
pling propagator t/ sll „ is applied. In order to preserve tr[8H 2 ] we take 
e/2 for each propagator. The dashed lines enclose a PAREC block. 
Repetition of PAREC between two gates means applying repeatedly 
one PAREC block after another. 



4 Accuracy improvement using PAREC 

In this Section we address the issue of quantum error correc- 
tion (QEC) of errors induced by static imperfections during the 
algorithm. The random errors in gates can be corrected up to 



a certain reasonable limit by usual QEC schemes which how- 
ever require a significant increase of the number of qubits [l] . 
Here we study a different scheme to correct the effects pro- 
duced by the static imperfections propagator U m . One possible 
way to correct errors produced by residual static couplings was 
introduced in 1131 . The idea is simple: contrary to random er- 
rors the static imperfections lead to a coherent accumulation of 
errors |[T2l . If some randomness is introduced then the effect 
of the residual couplings changes each time and does not ac- 
cumulate coherently. The PAREC method |[T3l profits from the 
freedom of choice of computational basis and uses this freedom 
by conveniently changing repeatedly and randomly the compu- 
tational basis along the computation. In order not to change the 
algorithm and the desired result a special care must be taken 
to properly compensate for the changes made renumbering the 
computational basis. 




Fig. 4. (Color online) Success probability of IQPEA which deter- 
mines the phase with a precision of up to 10 bits, in the presence of 
static imperfections, as a function of the imperfection strength £ 2 . The 
curves show the influence of PAREC and how probability is enhanced 
with the increase of the number of times PAREC is applied: (□) No 
PAREC; (■) 1 time; (o)5 times; (•) 10 times. No random imperfec- 
tions are present (£j = 0). Averaging is done over 2000 randomly cho- 
sen phases. 



The procedure is represented schematically in Fig. [3] To 
change the computational basis we first pick randomly from the 
set of Pauli matrices and identity matrix {X[,Yi,Zi,Ii} (where 
; = 1 , 2), and apply them to each qubit. We suppose that the 
time it takes to apply the Pauli operators is much shorter than 
any other time scales. We keep the information of this first 
choice, Xi,Z 2 in the Figure, and implement the suitably trans- 
formed gate (Zi ®X 2 )U\ (Z\ <S>X 2 ). After that to come back to 
the original basis, the operator (Z\ ®X 2 ) is applied again. The 
places where f/ SBt has acted are represented in Fig. [3]by thick 
vertical gray (red) lines. This procedure is repeated before and 
after each gate, but of course the key is that a new random se- 
quence of Pauli operators is drawn, in the figure (Y\ ®Z 2 ). So 
it is clear that between U\ and U 2 the imperfection propaga- 
tor U m has acted on different bases. We have called the com- 
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plete sequence of Pauli operators that act between two gates 
one PAREC block. 

The effect of PAREC can be seen in Fig. [4] The black solid 
curve (with □ symbols) shows the success probability when 
static (but not random) imperfections act. The gray (red) line 
(with ■ symbols) show the result when one PAREC block is 
applied between each gate of the algorithm that already gives a 
considerable gain. 

If instead of one PAREC block we introduce many of them 
(iVpAREc) keeping At fixed, and supposing that the time to imple- 
ment the Pauli gates is negligible, then the imperfection Hamil- 
tonian can be described with the help of the transformation 



V^VpAREC V2A^AREC 

As a consequence, the coherent effect of static imperfections is 
suppressed. This can be seen in Fig. [4] The gray (red) curves 
show the success probability as a function of £2 for different 
values of N PAREC (up to AWec = 10). As N PAREC grows the prob- 
ability grows accordingly. In the ideal limit of infinitely many 
PAREC blocks between gates (with a fixed gate-to-gate time) 
the success probability tends to constant maximum value for 
all e, a result which reminds us of a Zeno-like effect [ 16 1. This 
is also illustrated in Fig. [3] (top), where the dark region in the 
density plot of the success probability indicates the limiting 
value attained for large values of N ?AREC . The maximum value is 
the ideal value with perfect gates which is only limited by the 
value of the reminder 8 defined in ©. Nevertheless, the limit 
can be attained only theoretically because the time between IQ- 
PEA gates cannot be fixed if we add (ideally) infinitely many 
PAREC gates, no matter how fast we can implement them. 

Up to now we have considered PAREC with perfect Pauli 
gates while now we turn to a more realistic situation. With 
this aim we also consider the possibility of random phase er- 
rors of IQPEA gates to be also present in the PAREC Pauli 
gates. Therefore we expect that in the presence of random im- 
perfections, both in the IQPEA and in PAREC, there will be an 
optimal value of N PAREC after which the presence of too many 
faulty Pauli gates yields PAREC useless. This is demonstrated 
in Fig. [3] (bottom). For this a further consideration must be 
made. In Fig. [2] for illustration reasons only, we took E\ = £2 
such that the strength of both effects is comparable. However, 
we expect that in experiments the effect of random phase errors 
can be reduced to a minimum, so that in fact we have to assume 
£2 > £1 • For the plot in Fig.|5]we took £2 = 5 £1 . As a result the 
maximum of the success probability occurs approximately at 
^Vparec = 5 after which the PAREC method looses its efficiency. 
The position of the peak as well as its height depends on the ra- 
tio £ 2 /£i . The obtained data show that the IQPEA with PAREC 
can operate reliably even in presence of relatively strong static 
imperfections. 




Fig. 5. (Color online) Success probability (shown by color) of the IQ- 
PEA which determines the phase with a precision of up to 2~ 10 as a 
function of the number of times N PAREC PAREC is applied and of the 
static imperfections strength Top: the case where only static im- 
perfections are considered (£1 = 0). Bottom: both static and random 
imperfections are present (including in PAREC gates); here random 
errors in gates are also present, their strength is taken as E\ = £2/5. 
Averaging is done over 2000 randomly chosen phases. 



qubits. We have shown that static imperfections produce a dra- 
matic drop of success probability even for algorithms involving 
a rather small number of gates. In this context we have tested 
the PAREC method [ 1 3 1 and shown that it improves signifi- 
cantly the computation accuracy, even if the method is more 
suited to algorithms with a larger gate sequence involved. We 
also present results with repetitions of the PAREC method that 
produces a Zeno-like effect in preservation of probability. Even 
though, the realistic scenario would suggest a small N PAREC (may 
be even A^parec = 1), the results obtained demonstrate a convinc- 
ing improvement of the algorithm success probability induced 
by PAREC. The extension of the IQPEA circuits to a larger 
number of qubits is straightforward, as well as the PAREC im- 
plementation. 



5 Summary 

To summarize, we tested the effects of static imperfections in 
the IQPEA 1 6 1. Due to its simplicity this algorithm can be used 
as a benchmarking circuit for quantum computers with two 
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